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Abstract 

We show that higher-order coefficients required to perform threshold re- 
summation for electroweak annihilation processes, such as Drell-Yan or Higgs 
production via gluon fusion, can be computed using perturbative results de- 
rived in Deep Inelastic Scattering. As an example, we compute the three-loop 
coefficient D^^\ generating most of the fourth tower of threshold logarithms 
for the Drell-Yan cross section in the MS scheme, using the recent three-loop 
results for splitting functions and for the quark form factor, as well as a class 
of exponentiating two-loop contributions to the Drell-Yan process. 
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1 Introduction 



Soft gluon resummations El Ej have proven to be a valuable tool in perturbative 
QCD. They have provided a deep understanding of the structure of perturbation 
theory to all orders, which has in turn opened the door to studies on nonperturbative 
effects, and they have also been extensively used in phenomenology, broadening 
the range of QCD predictions towards the edges of phase space, where even hard 
processes are dominated by multiple soft gluon radiation. 

Resummation is closely related to factorization [4 . For threshold resummations, 
the hard partonic cross section for a given QCD process can be expressed as a con- 
volution (with respect to the energy fraction carried by hard partons, x) of different 
functions responsible for soft, collinear and hard radiation. The convolution turns 
into an ordinary product upon taking a Mellin transform. Logarithmic enhance- 
ments as X ^ 1 turn into logarithms of the Mellin variable N, and these logarithms 
can be shown to exponentiate, using evolution equations for the various functions 
involved in the factorization. 

To be precise, the resummed exponent is expressed in terms of moments of 
distributions singular as x ^ 1, 

as well as terms independent of A^, corresponding to moments of 6{1 — x) [3^. The 
pattern of exponentiation is nontrivial: in general, a perturbative calculation will 
contain terms of the form log^'^ multiplying the Born cross section, whereas 
in the exponent one finds at most terms of the form log'^'*'^ A^. Furthermore, a 
g-loop resummed calculation will determine completely the coefficients of the terms 
in the exponent proportional to log'^"'"^"^ A^, to all orders in a^. Such terms are 
usually described as N^~^LL, with leading logarithms (LL) determined at one loop, 
next-to-leading logarithms (NLL) determined at two loops, and so forth. 

Recently, the scope and expected precision of a range of QCD calculations have 
been extended in a remarkable series of papers by Moch, Vermaseren and Vogt 
(MVV), who computed first the three- loop contribution to the QCD splitting func- 
tions [nilZj, and then the complete three-loop DIS coefficient functions |8j, in what 
is arguably the most complex perturbative calculation ever carried out in quantum 
field theory. Their results both test and extend the range of threshold resummation 
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for DIS, which can now be performed exactly to N^LL accuracy. Furthermore, N^LL 
terms can also be determined, up to a single unknown coefficient requiring a four- 
loop calculation, the fourth-order contribution to the cusp anomalous dimension of 
a Wilson line in the MS scheme. It can, however, be argued convincingly that the 
numerical effect of this coefficient is negligible p . Thus soft resummation for DIS 
can now be tested at the level of the fourth tower of logarithms, providing nontrivial 
checks on the convergence of the expansion as the logarithmic accuracy is increased. 

Another class of benchmark cross sections for soft gluon resummation is given 
by electroweak annihilation processes in hadronic collisions, comprising Drell-Yan 
dimuon production, electroweak boson production, and Higgs production via gluon 
fusion. The inclusive cross sections for these processes are known to NNLO fOl El 
IT^ . and with the knowledge of the three- loop splitting functions the corresponding 
resummation can now be performed exactly at N^LL level, both in the MS and 
in the DIS factorization schemes. Lacking a three-loop calculation, however, N^LL 
terms are unknown, except for running coupling effects. It is the purpose of this 
letter to show that, using only results extracted from the three-loop DIS calculations 
of MVV, as well as known two-loop perturbative results for electroweak annihilation, 
one can bring the accuracy of threshold resummation for these processes in line with 
DIS, performing N^LL resummation up to the unknown, and very likely negligible, 
contribution of the four-loop cusp anomalous dimension. 

In the following, we will concentrate on the Drell-Yan cross section in the MS 
factorization scheme, although the reasoning is readily generalized to other elec- 
troweak annihilation processes and to the DIS scheme. We will make use of a 
factorization derived in jS], where the complete exponentiation of A^- independent 
terms was proven, to show that the coefficients of single-logarithmic contributions 
at g loops in the resummed exponent are completely determined by the knowledge 
of the (yf-loop nonsinglet splitting function, simple poles in the (yf-loop quark form 
factor, and A^-independent terms dX g — 1 loops in the Drell-Yan cross section. We 
will explicitly compute these coefficients at the three-loop level, and provide a gen- 
eral ansatz for their expression to all orders. These results will be useful in refining 
the theoretical prediction for processes of great interest at the LHC, such as Zq 
production and Higgs production via gluon fusion, by extending our knowledge of 
soft-gluon effects, and our control of the theoretical uncertainty due to uncalculated 
higher-order perturbative as well as nonperturbative corrections. 
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2 Factorization and exponentiation 



Our starting point is the unsubtracted partonic cross section for the Drell-Yan pro- 
cess. Near partonic threshold, its Melhn moments can be factorized as E] 

uiN,e) = |r(Q2,e)|2 iMN,e)f [/^(iV, e) + 0(l/iV) . (2.1) 

Here ipnlN, e) is the Melhn transform of a quark distribution at defined energy frac- 
tion, responsible for coUinear divergences, Uii{N, e) is an eikonal function describing 
the effects of soft gluon radiation at large angles, and r(Q^, e) is the (timelike) quark 
form factor. Near threshold, where all gluon radiation is soft, the quark distribution 
obeys a Sudakov-type evolution equation which can be solved in exponential form, 
as 

MN, e) = exp I £ dz^j'^ {a ((1 - y)'Q') , e) | . (2.2) 
Similarly, eikonal exponentiation applies to the soft function Ur, which can be writ- 

Ur{N, e) = exp I - dz ^ gu [a ((1 - zfQ^) , e) | . (2.3) 
The electromagnetic quark form factor F, on the other hand, is defined by 



r/.(pi,P2;/^ ,e) = (0|J^(0)|pi,p2) = -iecg v{p2)l^,u{pl) V[Q , (2.4) 

and it is one of the best understood amplitudes in perturbative QCD. Its logarithmic 
dependence on the scale can be determined using renormalization group and 
gauge invariance [ini El CSj) and the resulting evolution equation can be solved 
explicitly in dimensional regularization I6j, yielding the exponential expression 

^(Q^ e) = exp |i j^''^ ^ \yK (a., e) + G {a {e) {a (A^ 

(2.5) 

where 7E'(as) is the cusp anomalous dimension jHHIH], G(Q;s,e) collects all other 
scale-dependent terms, and is finite as e ^ 0, while K{as-, e) is a pure counterterm. 
A key feature of Eqs. ()2.2|l - ()2.5j) is the usage of the dimensional running coupling 
a(^^), defined in (i = 4 — 2e by the equation 

(3{e,a) = -2m + m , (3{a) = -—yhJ-) , (2.6) 

where 6o = (HCa — 2n/)/3 and hi = {17C\ — 5Ca^/ — 3Cf?2/)/6 in our normal- 
ization. Through a, integration over the scale of the coupling generates all infrared 
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and collinear poles in Eqs. ()2.2p - ()2.5|) . so that all functions appearing in the expo- 
nents are finite as e — 0, with the exception of the counterterm K in the quark 
form factor, whose only effect however is to cancel singularities arising from the 
^-independent limit of integration in the integral of the anomalous dimension 7^^. 
Further, dimensional continuation of the coupling regulates the Landau pole, which 
lies on the integration contour in c? = 4, allowing for an explicit evaluation of the 
exponents in terms of analytic functions of as and e ^01 ISHj • 

Our next task is to perform mass factorization on Eq. ^2.\\ . We do it here in 
the MS scheme, where we can make use of the expression jl] 



'MS = exp 



(27) 

Here A{as) can be extracted from the singular behavior of the nonsinglet QCD 
splitting functions as z — > 1, and is known to be related to the cusp anomalous 
dimension by A{as) = •jKicis)/'^, while Bs{as) is the coefficient of 5(1 — x) in the 
same splitting function. Once again, it is easy to see that (A^, e) is a pure 
counterterm, with all poles generated by integration over the running coupling. 
Clearly, Eq. ()2.7|1 is a simple exponentiation of the splitting function in the IR limit, 
including running coupling effects. Since it does not have an obvious diagrammatic 
interpretation (see, however, Ref. there is a certain amount of arbitariness in 

distinguishing real and virtual contributions in Eq. (|2.7|) . This arbitrariness was 
exploited in Ref. [S] to define 

</.j^(iV,e) = 0v(e)<^/?(iV,e), (2.8) 

where 



0y(e) = exp 




K (a., e) + G{a (e)) + \ /f '^Ik (« (A 



• (2.9) 



The structure of Eq. ()2.9|) clearly mimicks that of the quark form factor, Eq. ()2.5p . 
and in fact it is designed so that 0y(e) will precisely cancel all IR and collinear 
poles arising from T{Q'^,e). This requirement, together with the requirement that 
(/)y(e) be a pure counterterm, uniquely fixes the new function G(as), which can be 
determined recursively from G{as-, e), as was done explicitly in Ref. We are now 
ready to give our final expression for the Drell-Yan partonic cross section in the MS 
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scheme, which is 



u;{N,e) 



l^(g^6) 



(2.10) 

This expression has the important feature that virtual and real contributions are 
separately finite. Factoring out the virtual part (N) = |r(Q^, e)^ /{(j)y{e)y , and 
mapping the real terms to the conventional expression for the resummed Drell-Yan 
cross section in the MS scheme, including iV-independent terms as done in Ref. 
we are lead to our basic equation 



exp 



{MN,e)y 



[as) + 



dz ■ 



I- z 

+ D («, ((1 - z)2g2)) I +o{i/N). 



^(«.(/i^)) 



(2.11) 



Eq. (j2.1H) spells out our basic strategy to determine the resummation coefficients: 
(iV) must be finite by the factorization theorem, given our construction of the 



(R) 
MS 



virtual part Q^{N); the poles arising from the denominator, furthermore, are 
completely determined by the splitting functions and by the quark form factor, as 
seen from Eqs. (12. 7|) and ()2.9|) : requiring their cancellation determines a subset of the 
perturbative coefficients of the numerator functions, which are sufficient to control 
the expansion of the functions A and D. 



3 Constraints from finiteness 

The scale dependence of u)j^ (N) can be explicitly computed order by order making 
use of the exponential expressions for the functions ipR, Ur and (pR. An important 
point is the fact that ipR and Ur are renormalization group invariant P], which 
determines explicitly the scale dependence of their exponents. Consider for example 
the quark distribution ipR. Imposing RG invariance leads to 

+ ^) (^^^' = ' (3-1) 
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which can be solved perturbatively using the exphcit expression for the (3 function, 
Eq. ()2.6|) . and writing 



n=l 



«s, e) = II ( ^ ) «^Jr^ e 



(3.2) 



where from now on ^ will denote the ratio of the relevant scale (here (1 — x)Q) to 
the renormalization scale, for which we take fi = Q. Alternatively, one can impose 



(^, as, e) = (1, a(0, e) = 



n=l 



TT 



(3^3) 



which also determines the scale dependence of the perturbative coefficients k^'* (^, e). 
Using for the running coupling the solution of Eq. ()2.6|) expanded to three loops 



a as, e 



47re 



(3.4) 



one finds 



4'^(e,e) = 4'^(l,e)r^^ + ^.«(l,e)r^^(e--l) , 

= 4'^ (1, e) + Lf (1, e) + (1, e)) {C^' - l) 



(3.5) 



■^4^(M)(|-..)rne--i 



(3.6) 



with analogous results holding for the function gu{C,,as,e). The last formal step is 
to use the finiteness of and gu as e ^ to expand the e-dependent coefficients 
as 

00 00 
4Hl'^) = E4i^' ' 9lj\h^) = Y.9^U'^ (3-7) 



fc=0 



fc=0 



as well as 



G(a.,6) = EG(^He)(-)' =E(-J T.Gi 



(p) gfe 



(3.8) 



^ ^ p=0 ^ ' k=0 

Expanding, in a similar way, the various other functions involved in Eq. ()2.10p in 
powers of tts/vr, one can easily determine the structure of IR-collinear poles, by 
computing simple integrals. 
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It is instructive to briefly examine the information tliat can be extracted at tlie 
one-loop level. From Eq. ()2.1H) one derives 



,,(1) ^ ,.(1) 



9u,\ + ^,/,,2 \ 

2 J 

= Fig + D^^^Vq{N) + 4A(^)pi(Ar) . (3.9) 

The cancellation of double poles requires, unsurprisingly, that k^^q = 7^^''. Cancel- 
lation of single poles yields two equations, since the coefficient of the distribution 
Vo{N) must separately vanish. One finds that A^^'' = 4^o/2 = 'y^^ /2 (the factor of 
2 being a matter of historical conventions); further, one finds that a combination of 
coefficients of Ur and ipn is determined by yielding 

+ = -45?^ + 2G« . (3.10) 

Turning our attention to finite terms, we see first that the coefficient of the lead- 
ing distribution Vi{N) is confirmed to be A^^^ = 'y^j^ /2: had we not assumed the 
function A{as) appearing in to be the same as the one featuring in the resum- 
mation, this result would now have been derived at one loop. Next we see that 
single logarithms are given by the same combination of Drell-Yan coefficients that 
was determined by the cancellation of simple poles. This determines D^^^ in terms 
of DIS data as 

^(1) _ 4^(1) _ 2G(i) . (3.11) 

Finally, the one-loop exponentiated constants are given by = {gj)\ + 4\)/2. 

Clearly, all the coefficients involved at one loop are known or easily computed. 
For example one finds P], in the MS scheme. 



(1, 6) = 2Cpe^^- ~ '\ , (1, 6) = -2C^e^^- ~ '\ 

V> V , ; r(2-2e) ^ ^ 1(2 - 2e) 



while, as derived in [F, G^^) = c'i^ = 3Cf/2. It is well-known that sj^^ = 3C^/4, 
so one finds consistently 

^^^^ = , F« =-^C(2)C^, (3.13) 
as confirmed by a direct one-loop calculation of the Drell-Yan cross section. 
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At two loops, the pattern repeats iself with a few twists. The cancellation of triple 
and double poles brings in no new information, except the fact that the function 
begins to differ from by running coupling effects, 

^(2) _ ^(2) , ^0 / (1) 3 (1) \ (2) 1 , ^ 1 

'^V,o - 7x + y \ 9ufi + 2^^,i J -^K + 2"o<-^F • 

This however is just a reshuffling between ipji and Un, in fact at the level of single 
poles the effect cancels and one finds, as expected, that requiring the cancellation 
of VQ{N)/e terms yields A*^^) = '~fK^/'2 [221 123 • A^-independent single-pole terms, on 
the other hand, constrain a combination of coefficients of gu and k^, namely 

(2) 

+ ^ = -4Bf ' + 2G<^I + ^ {aii\ + |.»> ) . (3^15) 

Turning to finite terms, one finds that once again running coupling effects involving 
tpR and Ur cancel, and single logarithms are determined by 

D(^) = 45f ) - 2G(^) - ^ {gS\ + ) = 45f ^ - 2&^^ - . (3.16) 

All required ingredients are known: from Refs. PIES], while G^^) = G^q^ - 

BqG^i^ /A was given in One finds then 

D^'^ = (-^ + yC(2) + ^C(3)) CaCj. + - ^C(2)) njCj. , (3.17) 

which agrees with a direct comparison jH 12^] with the two-loop calculation of 
Ref. [TU], in the spirit of [2Z|- Exponentiated two-loop constants are also constrained 
by 

p(2) _ 1 I A2) i^^A bo / (1) 3 
where running coupling effects are readily evaluated using Eq. ()3.12|1 . 

4 The coefficients D^'^^ at higher orders 

It is straightforward to continue the analysis at three loops. As expected, the cancel- 
lation of quartic and triple poles at three loops in Eq. 1)2.111) is achieved automatically 

^Notice however a misprint in Eq. (4.6) of Ref. the coefficient of CaCf in Gq^'' should read 
(2545/108 + llC(2)/3 - f 3C(3)))/4. 
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as a consequence of lower-order constraints. Double poles specify the relationship 
between and at the three-loop level; using Eq. ()3.15|) one can write 

^(3) _ (3) ^ (2) _ &| (1) , / (1) 3 

As before, running coupling effects do not affect the known relationship between 
A{as) and 7i^(as): demanding the cancellation of of Vo{N)/e terms at this order in 
fact yields A^^^ = j^^^ /2. A^-independent single-pole terms, on the other hand, yield 
the constraint 

(3) 

sSl+"-f = -4Bf + 2G«) + ^(jgl + 54;)) 

The finite coefficients of Vi{N) with i = 1, 2, 3 provide nontrivial tests of the results 
achieved so far. Further, concentrating on single logarithms, and using Eq. (j4.2j) . 
one finds that 

= 4Bf-2G^')-i,„fi-|F« . (4.3) 

The detailed structure of the coefficients in terms of the functions gu and k^, as 
before, turns out to be irrelevant, and the aswer is simply expressed in terms of lower 
order contributions to the function Fjj^ {as)- This is remarkable, but easily under- 
stood: in fact the details of the factorization given in Eq. ()2.H1 . while conceptually 
crucial to prove formally the exponentiation of logarithms to all orders, cannot affect 
the overall structure of IR-collinear poles: one could, for example, define a modified 
quark density including eikonal effects, and poles would still cancel. Inspection of 
Eqs. ()3.11|) . ()3.1(i|) and ()4.3j) leads us then to the following all-order ansatz for the 
function D[as), which summarizes the results of our work. 

D{as) = ABsias) - 2 + -^F^ M ■ (4.4) 

The function D{as), governing threshold resummation for electroweak annihilation 
at the single-logarithmic level, is thus completely determined at order a" by the 
knowledge of virtual contributions to the nonsinglet splitting function, and IR- 
collinear poles of the quark form factor, to the same order, plus the value of ex- 
ponentiated A^-independent terms arising from real emission at order 
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We are now in a position to evaluate the three-loop contribution to the function 
D{as), thanks to the recent results of MVV. The three-loop contribution to the 
function Bs{as), in fact, is given in Ref. the three-loop coefficient of the function 
G{as) is given (in [5,) by the expression 



= - ^Gf ) - + , (4.5) 

4 4 lo 

and all relevant coefficients in the expansion of the function G{as, e) can be found 
in Ref. |2H1, where MVV use their results for DIS structure functions to evaluate 
explicitly the quark form factor at three loops; finally, the value of Fy^ (a^) at two 
loops can be extracted by comparing our exponentiated expression with the two-loop 
calculation of Ref. jTU]. We find 

^ - (il->)nc^(^)->))c.c. 

Collecting all ingredients, or result for D^^^ is 

V 23328 324 ^^ ^ 60 ^ ^ ^ 108^^ ' 6^^ ' ^\ ) ) a t 

/31313 1837^,^, 23^2/oN 155^, ^ ^ 

+ C(2) + — C (2) C(3) UfCACp 

V11664 324 ^^ ' 30^ ^ ^ 36 V ^ 

+ ( w - > 4^^<^) - >') "^^^ 

The coefficient of the highest power of in D^^^ can be independently checked by 
comparing it with the renormalon calculations of [23 and j^: indeed, their results 
agree with the last line of Eq. ()4.7|) ^. 



5 Discussion 

We have analyzed threshold resummation for the Drell-Yan process in the MS 
scheme, in light of the recent results obtained for Deep Inelastic Scattering by MVV. 
Building upon a factorization proposed in Ref. 0, we have been able to derive a 
■^We thank Einan Gardi for pointing out this check to us and providing us with his results. 
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general relationship expressing the function D{as), responsible for threshold loga- 
rithms in the Drell-Yan cross section at single-logarithmic level, in terms of data 
requiring the knowledge of the virtual part of the nonsinglet splitting function, and 
the singular terms in the quark form factor, at the same perturbative order, plus 
a well-defined set of A^-independent terms arising in the Drell-Yan cross section at 
lower orders. Our main result is Eq. ()4.4|) . and, using MVV results, it has enabled 
us to evaluate the three-loop coefficient D^^\ given in Eq. ()4.7|1 . 

An immediate question is whether our results extend to the case in which the hard 
annihilating partons are gluons, which is relevant for the process of Higgs production 
via gluon fusion, in the effective theory with the top quark integrated out. It is, 
in fact, easy to show that an equation identical in form to Eq. ()4.4|) holds also for 
gluon-initiated electroweak annihilation, provided the various functions involved are 
appropriately redefined: in fact, threshold resummation in that case can still be cast 
in the form of Eq. ()2.1H) . with 2A{as) replaced by the cusp anomalous dimension for 
a Wilson line in the adjoint representation, 274^(0;^), and two new functions Dg{as) 
and F-^ (a^). The MS distribution can be similarly defined for initial gluons, 
with Bs{as) replaced by the virtual part of the appropriate gluon splitting function. 
The gluon form factor obeys an exponentiation identical in form to Eq. ()2.5|) . All 
ingredients are thus in place to yield Eq. ()4.4|) . A more delicate question is whether 
this implies a simple relationship between the perturbative coefficients of D and Dg. 
Up to two loops, one verifies by explicit calculation jT^I that Dg can be obtained 
from D by simply replacing the overall factor of with Ca, just as one does in 
deriving Ag from A. It is unlikely, however, that such a simple behavior will persist 
to all orders: in fact, while it is natural to expect that purely eikonal quantities such 
as A or the function Ur will be sensitive only to the representation of the gauge 
group in which the eikonal line is placed, not all information encoded in Eq. ()4.4|) 
arises from eikonal lines; it is known, for example that sub leading poles in 
the gluon form factor cannot be obtained from the quark form factor with such a 
simple prescription. Even eikonal functions would probably require a more careful 
treatment at high enough order, when high-rank Casimir operators constructed out 
of the symmetric SU (N) tensors dabc come into play. 

All this notwithstanding, we argue that at the three-loop level the simple pre- 
scription is still valid, and one can in fact compute D^^^ by simply replacing the 
overall factor oi Cp with Ca- To see it, one can make use of an observation of 
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Ref . jSI] , already exploited in Ref . . According to this observation, it is possible 
to isolate in the quark form factor, and specifically in the function G{as, e), a class 
of maximally nonabelian contributions, dubbed fj^"^'^^ in Ref. [221, which exhibit the 
same behavior as the eikonal anomalous dimenson A {i.e. they obey the simple re- 
placement rule, as verified up to three loops in |32 ). We have explicitly checked up 
to three loops that in fact the leading terms of our equation, ABg'^ — 2G'^^\ coincide 
with the maximally nonabelian factors fl up to an irrelevant multiplicative factor. 
Since the remaining term in our Eq. ()4.4j) is a running coupling effect, determined at 
lower orders where the replacement rule is known to apply, we conclude that indeed 
D'^g^ is also given by Eq. ()4.7|) . with the overall Cp replaced by Ca- 

We conclude by noting that we expect these results to be useful for hadron 
collider phenomenology. In fact, along the lines of j^, the knowledge of /^'■'^^ allows 
to perform N^LL threshold resummation for Drell-Yan and Higgs production, to 
what is expected to be a very good approximation. This can be used not only to 
provide a more accurate QCD prediction for these processes, but also to check for 
the stability and the convergence properties of both ordinary perturbation theory 
and the expansion of its resummed counterpart in towers of logarithms. Finally, we 
note that several of the building blocks of our analysis also enter in resummations 
and high-order perturbative calculations for more complicated processes at hadron 
colliders (see for example [22] )• It would be interesting to study the extent to which 
our techniques can be apphed also in that context. 
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Note added 

While our paper was being written, S. Moch and A. Vogt completed their own cal- 
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culation of D^^\ for both quark- and gluon-initiated scattering jSlj, using a different 
line of argument. Their results completely agree with ours. 
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